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We have calculated high temperature expansions for the momentum distribu- 
^ ■ tion function n k and the equal time spin and density correlation functions <S(q) 



and iV(q) of the two-dimensional t-J model. On extrapolation to low tem- 
peratures we find that n k has a step-like feature at k F and <5(q) has 2k F as a 
characteristic wavevector, whereas iV(q) has 2k F F as a characteristic wavevec- 
tor. Here k F and k F F are the Fermi wavevectors of the nearest-neighbor square 
lattice tight-binding and spinless fermion models, respectively. By comparison 
to the known results for one dimension this suggests spin-charge separation 

in the two-dimensional t-J model. 
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The study of strongly correlated electrons in two- dimensions (2D) is currently one of 
the most interesting and controversial topics in condensed matter physics, particularly with 
regard to high temperature superconductivity (HTSC) HJ. Anderson has put forward 
the idea that the ground state of strongly correlated electron systems in 2D is a Luttinger 
liquid analogous to the case in one dimension. In ID, a Luttinger liquid has spin and charge 
degrees of freedom with different velocities and wavevectors, behaving at low energies as 
independent elementary excitations, a situation which has become known as spin-charge 
separation @. 

Determining whether or not spin-charge separation can also occur in 2D has proven 
quite difficult. We have calculated high temperature expansions for equal time correlation 
functions (ETCF) of the 2D t-J model || to investigate this possibility. We find two distinct 
characteristic wavevectors for the spin and charge degrees of freedom, 2k F and 2k^ F defined 
below. This shows that the spin and charge degrees of freedom have different distributions 
in the Brillouin zone and provides evidence for spin-charge separation in this model. 

We consider the 2D t-J model on a square lattice, where the Hamiltonian is 

H t J = "* E iAaCja + h.C.) + J E S, " S„ (1) 

<ij>,a <ij> 

with the constraint of no double occupancy. The constraint represents the strong correlations 
between the electrons and is difficult to treat by conventional many-body techniques. 

We have studied three ETCF of this model using the high temperature series expansion 
method. These are the single spin momentum distribution function, n k , and the equal time 
spin and density correlation functions, 5*(q) and iV(q), defined by the relations 

n k — E e * k r ( c Qo- c ro-)' 

r 

S(q)=Ee iq - r (^>, 

r 

iV(q)=E eiq ' r ( A ^oAn r ), (2) 

r 

where the angular brackets refer to thermal averaging in the grand canonical ensemble, 
$r = \ TjolP c \ a a a3 c ri3 an d An r = J2a c L- c ro- — Ti. Here u is the average density of electrons. 
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The expansions are calculated for n k ||] to eighth order and S(q) [|J and N(q), to tenth order 
in the reciprocal temperature T -1 , first at fixed fugacity, and then by a change of variables 
at fixed n. The series are extrapolated at fixed n by Pade approximants to determine the 
low T properties. 

The form of the ETCF for tight-binding (TB) (non-interacting, spin-half electrons with 
nearest neighbor hopping on a square lattice) and spinless fermions (SF) (physically, SF 
are fully spin-polarized TB electrons, freezing out the spin degrees of freedom and doubling 
the number of occupied k states) in 2D are helpful in understanding the t-J model results 
presented below. They are given by 



where g = 2 for TB or g = 1 for SF, and 4S'(q) = iV(q) for TB. From the form of this 
equation we can see that at T = and n < 1 the ETCF will saturate at n when n k and 
n k+q no longer overlap (note that for SF with n > 0.5, iV(q) saturates at 1 — n when the 



existence of a sharp Fermi surface. 

In Fig. ^ we compare iV(q) of the 2D t-J model to iV(q) of SF at the same density for 
T/J = 0.5, J/t = 0.5 and a range of n. The similarities are remarkable throughout the 
Brillouin zone, with the differences near the T point due to the t-J model having a larger 
compressibility than SF. To focus the discussion below we now limit ourselves to two sets 
of parameters outside of the phase separated [JJ or ferromagnetic || regions of the 2D t-J 
model. We fix J/t = 0.5 and consider n = 0.75 and n = 0.20. The results for n k , S(q) 
and iV(q) along q VM = (0, 0) — > (tt, it) are shown in Fig. |. We see that n k ~ 1/2 at k F rM , 
the Fermi momentum of the TB model at the same density [|]|| and S(q) is enhanced over 
its TB value and either flattens out or has a peak at q w 2k F rM . However, the most 
anomalous curves are for iV(q). They are suppressed from their TB values and flatten out 
at q ~ 2k" F F rM , the Fermi momentum of SF at the same density. We observe no feature at 
q = 2k F , though N(q) flattens out more gradually for n = 0.20 than for n = 0.75. 




(3) 



hole Fermi surfaces no longer overlap). The Kohn anomaly |J at 2k F or 2k F F is due to the 
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For comparison, we recall the behaviors of n k , S(q) and N(q) for the U/t — > oo Hubbard 
model (J/t — > t-J model) in ID. From the work of Ogata and Shiba |K| we know that n k 
has a power law singularity at k F with nk F = 1/2. Also S(q) has a peak at 2k F and N(q) 
has 2k s F F = 4k F as a characteristic wavevector. For arbitrary U/t, N(q) also has a feature at 
2k F due to a mixture of spin and charge excitations, but the 2k F F feature is due to charge 
alone and shows that the charge degrees of freedom truly reside at k s F F . In addition to 



the singularity at k F , n k has a singularity at 3k F , but with a very small step size. 

Our calculated 2D ETCF show behaviors very similar to their counterparts in ID. The 
characteristic wavevectors for S(q) and iV(q) are 2k F and 2k F F , respectively, which we 
believe implies low energy spin degrees of freedom near k F and low energy charge degrees of 
freedom near k F F . Note that in 2D k F and k F F are incommensurate wavevectors; the charge 
degrees of freedom do not occur at a harmonic of k F , but at an independent wavevector. 
In Fig. [| we show k F and k F F for the whole Brillouin zone at n = 0.75 with representative 



nesting wavevectors. For weak coupling calculations of the 2D Hubbard model [I2[ and 
Gutzwiller projected free electrons [T3|] the picture is quite different. In these cases while 
S(q) is enhanced and iV(q) is supressed, they both have 2k F as a characteristic wavevector 
which is not what we find for the 2D t-J model. The behavior of n k is also similar to ID. 
The step in n k at n = 0.20 is comparable in size and shape to the TB model at the same T, 
but at n = 0.75 the step is much weaker and too smeared out to be explained by thermal 
broadening alone |4j]. We have not seen any evidence for a singularity at 3k F in 2D. This 
could be due to the relatively high temperature T / J = 1.0 in our calculation or possibly the 
angular averaging in 2D which is not present in ID. 

In ID the statistics of the excitations play no role, but in 2D they are important. Our 
data give no direct evidence on the statistics of the excitations in 2D, but we can formulate 
a hypothesis as to what they might be . If we think of a single electron as being composed 
of an elementary spin degree of freedom and an elementary charge degree of freedom, we 
would expect one of them to be fermionic and the other bosonic to give a fermionic electron 



14j| . Since n k shows a step at k F and S(q) has 2k F as a characteristic wavevector we assign 



the spin degrees of freedom as fermionic and the charge degrees of freedom as bosonic, but 
note that the charge degrees of freedom are not free bosons, but hard core bosons (HCB) 
to enforce the constraint of no double occupancy. This can be seen in Fig. [I] for n = 0.5 
where the data points near (ir, ir) are already more rounded than SF at T j J = 0.5. Further 
evidence for this point of view can be gained from the work of Long and Zotos [13] and 



Sorella, Parola and Tosatti [16 



We have also estimated the behavior of the HCB iV(q) by a flux phase mean field 
calculation. In 2D, HCB on a square lattice with nearest neighbor hopping can be exactly 
mapped into SF by attaching a quantum of magnetic flux 0o to each particle |l7j . If density 



fluctuations are not large we can replace the attached flux tubes by a uniform magnetic field, 
B = n0o, which will couple to the orbital motion of the particles. This corresponds to a 
SF model with a site dependent phase (the uniform flux phase). Using this flux phase mean 
field approximation we have calculated iV(q), with the results at T = shown in Fig. 2. 
The global features show a rounded flattening out of iV(q) at 2k^, F and general agreement 
with SF and the t-J model. For small q the approximation we are using gives iV(q) oc q 2 , 
but by general hydrodynamic arguments we know that for T = if the system has a finite, 
non-zero compressibility the q — > limit should be linear. The quadratic dependence is 
due to the "Fermi energy" of the flux phase sitting in an energy gap []18| |. Therefore the 



g 2 -dependence is an artifact due to our mean field approximation and should become linear 
after including fluctuations, which we will discuss in a future paper. 

More information on the interactions between the spin and charge degrees of freedom 
could be obtained by considering the 2D t-J model with a non-zero spin polarization. If 
the spin and charge are coupled we would expect both S(q) and iV(q) to change. However, 
if the spin and charge degrees of freedom are truly separate the characteristic q-vector of 
iV(q) should not be affected by a non-zero spin polarization but S(q) would now have 
transverse and longitudinal components with features at wavevectors that depend on the 
number of up and down spins. This has been observed by Ogata, Sugiyama and Shiba 
for the ID U — > oo Hubbard model. 



Having elementary degrees of freedom at k F and k|, F should have experimental con- 
sequences for the copper oxide planes in HTSC Neutron scattering experiments on 
La2_ x Sr x Cu04 show four incommensurate peaks centered around (n, n) that move with dop- 
ing. This can be understood in terms of the nesting properties of the weak coupling Fermi 
surface (521 where our results also put the spin degrees of freedom. The energy integrated 



weight of angle resolved photoemission is a direct measure of n k and photoemission has also 
been interpreted as supporting a large Fermi surface [fffi . But the transport measurements 



24fl are not so easy to understand in this picture. 

Reconciling experiments which show a large electron-like contour in k-space with a den- 
sity of n carriers, with transport measurements, which show a much smaller hole density of 
1 — 72 carriers, is one of the most puzzling problems of the copper oxides p4 |. Our results for 



the 2D t-J model show one way this might occur [@,^IJ. One expects the transport experi 



ments to couple most strongly to charge. For 1 - n < 1, the charge degrees of freedom at 
k|, F have a small hole-like locus in momentum space centered around (jr, tt) as shown in Fig. 
[3J. With this the transport measurements are satisfied. At the same time the spin degrees 
of freedom give a large, weak coupling Fermi surface also shown in Fig. [3J, which is seen in 
neutron scattering and photoemission experiments. We wish to emphasize that experiments 
which could probe iV(q) directly may prove to be very interesting for HTSC materials. 

In conclusion, we have studied the equal time correlation functions of the 2D t-J model 
by high temperature series expansion methods. We find that the spin and charge ETCF 
exhibit signatures of two different wavevectors: the characteristic wavevector for the spins 
being k F and that for charge k^, F , the Fermi wavevectors for TB and SF respectively. In 
comparison with the results for ID this suggests spin-charge separation in this strongly 
correlated 2D model. 
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FIGURES 

FIG. 1. Plot of N(q) at T/J = 0.5 and J/t = 0.5 along the irreducible wedge for a range 
of n. The data points are the t-J model and the solid lines are spinless fermions for the same 
temperature. The small vertical arrows are the T = locations of nesting vectors for spinless 
fermions. 

FIG. 2. Plots along the diagonal r — > M at n = 0.20 (a) Single spin momentum distribution 
function. Data points: t-J model, solid line: tight-binding model at T/J = 1.0; (b) Spin correlation 
function. Data points: t-J model, dashed line: T = tight-binding model; (c) Density correlation 
function. Data points: t-J model, solid line: T = flux phase mean field approximation for hard 
core bosons, dashed line: T = spinless fermions and dotted line: T = tight-binding model. The 
vertical dashed lines indicate the important wavevectors along this line in the Brillouin zone for 
tight-binding electrons and spinless fermions: nesting wavevectors 2k F TM and 2kf F rM , or Fermi 
wavevectors k F rM and k F F rM . (d) - (f) same as (a) - (c) with n = 0.75. K TM = (27r,27r) is a 
reciprocal lattice vector. 

FIG. 3. Fermi wavevectors for n = 0.75. Solid curve: tight-binding electrons, dashed curve: 
spinless fermions. The arrows are representative nesting wavevectors along TM and MX. 



11 



